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1.0 INTRODUCTIOK 


This document contains requirements for the asoent/desoent (A/D) navigation and 
attitude-dependent ancillary data products to be generated for the Space Shuttle 
Orbiter by the Mathematical Physics Branch/Mission Planning and Analysis Divi- 
sion (MPB/MPAD) in support of the orbital flight test (OFT) flight test re- 
quirements (FTR's), MPAD guidance and navigation (G4N) performance assessment, 
and the mission evaluation team. This document is intended to serve as the sole 
requirements control Instrument between MPB/MPAD and the A/D ancillary data 
users. The requirements preseti-ed herein are primarily functional in nature, 
but some detail level requirements are also included. 


2.0 ASCENT /DESCENT ANCILLARY DATA SUPPORT SCOPE 

The A/D ancillary data support for OFT mission evaluation activities shall be 
confined to providing postflight position, velocity, attitude, and associated 
navigation and attitude derived parameters for the Orbiter over the flight 
phases and time intervals shown in figure 1 . No ancillary data support related 
to the external tank (ET) or the solid rocket boosters (SKB's) shall be pro- 
vided. In addition, the A/D ancillary data products delivered shall be confined 
to those described in section 5.0 of this document. 


3.0 DATA PROCESSING APPROACH (FOR INFORMATION) 


This section summarizes the A/D ancillary data processing approach in terms of 
the functional flow diagram shown in figure 2. The approach is presented primar- 
ily for information purposes to acquaint prospective ancillary data users with 
the current data reduction plan. Elements of the approach are subject to change, 
pending further analyses of data processing techniques and procedures. 


3.1 INPUT DATA SOURCES AND INITIAL DATA PRODUCTS (TIER 1) 

Raw C- and S-band tracking data shall be obtained from the OPT Mission Control 
Center (MCC) through the facilities of the Ground Data Systems Division (GDSD). 
Selected onboard navigation and attitude telemetry data shall be obtained from 
the Institutional Data Systems Division/Data Reduction Complex (IDSD/DRC). 

Phototheodolite data shall not be processed in this task. 

3.2 INTERMEDIATE DATA PRODUCTS (TIER 2) 

Preprocessing or conditioning of the input tracking and telemetry data shall be 
accomplished by using interactive computer graphics technicjues . Ground and 
onboard clock times shall be tabulated to implement correlation of sensor data 
from the two input sources. Special event times will be provided by the JSC In- 
strumentation Integration Branch (WC6). 
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3.3 FINAL PROCESSING AND DELIVERED OUTPUT DATA PRODUCTS (TIERS 3, <4, AND 5) 

At the tler-3 level, an ephemeris solution (position, velooity, and accelera<- 
tlon) shall be generated first, using a navigation filter procedure that com- 
bines onboard and ground tracking data in the A/D flight regions area. (Both 
sources are available.) In regions where ground tracking is precluded (e.g,, 
during entry blackout), the solution shall be obtained by using onboard inertial 
measurement unit (IMU) data to propagate the state vector computed at the 
terminus of the navigation filter solution. Attitude data are then refined and 
time-correlated with the overall ephemeris solution. This result shall provide 
the basis for the delivered output products indicated in tier 4. 

Tier-4 processing shall Include the computation of the remaining navigation and 
attitude-related parameters needed to complete the parameter set established for 
the output product; i.e., the product intended to satisfy the ancillary data 
users. Ascent quick -look output data will be provided in this tier. The basic 
computation for the ascent quick-look outputs is the same as the final outputs 
with the exception of enhanced BET accuracy. The final descent best estimated 
trajectory (BET) is developed in tier 4 and reflects the meteorological data 
supplied to MPAD by WC6. 

MSFC via WC6 will provide the meteorological data for the ascent segment of the 
OFT flights. These data will be provided by MSFC 5 days after flight phase ter- 
mination. The meteorological data for the descent phase will be supplied by the 
National Meteorological Center (NMC) 3 weeks after the OFT flight via WC6. The 
NMC requires plots and tabulation of the descent ground track based on telemetry 
(TM) vectors for generating the meteorological data. A time history of center- 
of -gravity data for the ascent and descent segment of the OFT flights will also 
be provided by WC6. This information will be reflected in the ascent final and 
quick -look product. 

In association with the ascent and descent BET’s, a set (table 1(c)) of naviga- 
tion parameter accuracies will be generated as a special product that gives an 
estimate of the error associated with the tier-3 outputs, which includes the un- 
certainty associated with the observed meteorological data provided by the MSFC 
and NMC via WC6. 
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4.0 REPUIRED INPUTS 

Raw C- and S-band tracking data shall be obtained from the OFT MCC through the 
facilities of the GDSD. Selected onboard navigation and attitude telemetry data 
shall be obtained from the IDSD/DRC. 

MSFC via WC6 will provide the meteorological data and associated accuracies for 
the ascent segment of the OFT flights. These data will be obtained by MSFC 
through KSC. 

Meteorological data and associated accuracies for the descent phase will be 
supplied by the NMC 3 weeks after the OFT flights. The NMC requires plots and 
tabulation of the descent ground track (latitude and longitude) based on TM 
vectors for generating the meteorological data. Special event times will be pro- 
vided by WC6. 

In order that the task provide velocity and acceleration vector data with re- 
spect to the vehicle center of mass, the task requires as input a table that de- 
fines the location of the center of mass with respect to the navigation base. 

The definition of the center of mass with respect to the IMU navigation base 
will be coordinated through WC6. 


5.0 GENERAL ORBITER ASCENT /PESO ENT OUTPUT PRODOCTS 


Functional as well as some detail level characteristics are presented here for 
the Orbiter as cent/ descent ancillary output products. Both general and special 
ancillary output products are addressed in terras of the following output require- 
ment items: 

a . Parameters 

b. Coordinate systems 

c . Units 

d. Accuracies 

e. Time correlation 

f. Output data frequency 

g. Output data forms 

h. Data product delivery schedules 

i. Data product distribution 
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5.1 GENERAL OUTPUT PRODUCTS 

The output product described below is designed to provide Orbiter asoent/descent 
ancillary data to satisfy users with a single product. The product output param- 
eters and format will be common to all Orbiter ascent and, descent flight phases 
(tables I and II). 


5.1.1 Parameters, Coordinate Systems, and Units 

Table 1(a) provides the BET output product tape format, Table 1(b) describes 
the output tape header record, The product output parameters and their asso- 
ociated coor'^inate systems and units are given in table I(o). Definitions of 
the relevant coordinate systems are given in the appendix, together with equa- 
tions for selected output parameters to supplement their tabular descriptions. 
The product ephemeris parameters in table 1(c) all describe the motion of the 
Orbiter in navigation base and in terms of the center of mass, Table 1(d) de- 
scribes the BEr products listing format. 

There are some parameters that are output in several coordinate systems and have 
identical descriptions. These parameters are vector quantities whose components 
depend on the reference coordinate system. The following gives the generic de- 
scription and a more complete definition; 

a. Contact acceleration - The acceleration of the vehicle, excluding 
gravitational acceleration. 

b. Gravitational acceleration - The acceleration of the vehicle caused by the 
Earth’s gravitational force acting on the vehicle. The gravitational 
acceleration is computed by using a gravitational field model 
(appendix, p, A-21). 

c. Total acceleration - A vector sum of the contact and gravitational accelera- 
tions . 

d. Earth -relative velocity - The velocity of the vehicle with respect to the 
Earth-fixed coordinates (inertial velocity minus W x R, the Earth 
rotation vector crossed into the inertial position vector). 

e. Wind-relative velocity - The velocity of the vehicle with respect to the air 
mass through which it is moving. The vector is a result of subtracting the 
wind vector from the Earth-relative velocity vector. The wind vector, if 
not measured directly, is set to zero. 

Tables 11(a) through 11(e) presents the "one-time-only” parameters for block 
outputs. Table 11(a) presents the navigation block output format. The tape 
header record is given in table 11(b). The output tape parameter location is 
presented in table 11(c) with their respective description in table 11(d). 

Table 11(e) contains the output listing format. These parameters will remain 
constant during the duration of the mission phase. 


4 


78FM40 


5.1.2 Aoouracias 

The product generated by this task will satisfy the PTR accuracy requireraento to 
the extent possible baaed on the quality of the onboard and ground navigation 
sensor data available, 


5.1.3 Time Correlation 


Time correlation of the product output variables shall be accomplished as in- 
dicated in sections 3.2 and 3.3. 


5.1.*1 Output Data Freauenov 

Product output data shall be generated at one sample per second on the even sec- 
onds Greenwich mean time (GMT) and at requested event times provided by WC6. 


5.1.5 Output Data Forms 

Product output data shall oe ckslivered on Uni vac 1108 and CDC Cyber computer- 
compatible tape(s) (CCT's) accompanied by computer printout and microfilm 
listings. Tables 1(a) through 1(d) present the BET output product tape format, 
(table 1(a)), the header record (table 1(b)), the parameters, and associated 
word number and symbol (table 1(c)) as well as the computer printout listing 
format (table 1(d)). Tables 2(a) through 2(c) provide similar information as 
tables 1(a) through 1(d) but for the navigation block output parameters. 


5.1.6 Data Product Delivery Schedules 

The ascent quick-look product output shall be made available for distribution 
nominally 1 week after the flight phase. The final output product for ascent 
will be delivered 4 weeks after the flight phase and for descent, 5 weeks after 
the flight phase. All deliverables are dependent on receiving the required 
inputs. 


5.1.7 Data Product Distribution 

Distribution of the product output shall be accomplished as follows: 

a. WC6 shall serve as the sole distribution center for all product output. 

b. Requests for the data shall be made to WC6 where the appropriate JSC 
facilities will be directed to generate the required number of tape copies 
and/or listings of the tape(s). 

c. MPB shall provide WC6 with the appropriate tape reel numbers, a computer 
printout listing of the tape(s), and a microfilm of the listing. 
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d. MPB ahall provide IDSD with a copy of the aotual tape(s), 

e. MPB ahall maintain the reaponaibility for proper archiving of the 
aacent/deaoent ancillary data. 


5.2 SPECIAL OUTPUT PRODUCTS 

Special output produota are defined to be thoae requeated by a amall number of 
poatflight analyata whoae ancillary data requirementa are not aatiafied by the 
general output product. 


5.2.1 Parametera, Coordinate Syatema, and Unite 

Time hiatoriea of the estimated accuracy associated with a subset of the general 
product parameters of table I(o) shall be provided. The subset is identified in 
table 1(c), The accuracy estimates of the BET will reflect the uncertainties 
associated with the measured winds and atmospheric parameters. Plots and tabula- 
tions of the descent ground track will be provided to NMC via WC6 for the devel- 
opment of the entry wind and atmosphere profile. Another support item that will 
be provided under special output is a quick-look time history of the altitude 
and Earth-fixed velocity magnitude of the descent trajectory based on downlinked 
telemetry data. 

Coordinate systems and units of the latitude, longitude, altitude, and Earth- 
fixed velocity magnitude and navigation parameter accuracy outputs shall be iden- 
tical to those defined for the general product parameters. 


5.2.2 Ancurao ies 


The product generated by this task will satisfy the FTR accuracy requirements to 
the extent possible based on the quality of the onboard and ground navigation 
sensor data available. 


5.2.3 Time Correlation 

There is no requirement for the special outputs to be time correlated with re- 
spect to each other, but they shall be time tagged. 


5.2.4 Output Data Frequency 

The parameter accuracy estimates output rates will be at least approximately 
every 20 seconds. Outputs for the ground track and stripped telemetry parame- 
ters will be at the rate of the downlink data. 
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5.2.5 Output Data Forma 

Output data fopna will oonaiat of computer printout CCT plots based on the appro- 
priate user requirements, 


5,2.6 Data Product Delivery Schedules 

The ground track and stripped telemetry parameters will b® available 1 week 
after the flight. The accuracy estimates shall be distributed 3 w^eks after de- 
livery of the descent general product. 


5,2.7 Data Product Distribution 

Distribution of the discrete event special product shall be accomplished as 
discussed in section 5.1.7, 
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TABLE I.- ORBITER PRODUCT OUTPUT PARAMETERS 
(a) BET output products tape format 




Version 



CDC 

Univac 

Record 

Description 

1600 BPI 9 track 

800 BPI 9 track 

number 

of contents 

written by CDC 

written by CDC 

1 

Header record 804 

80 CDC S.P. word 

134 Univac S.P. 


characters of 

(60 bit) binary 

word (36 bit) 


Hollerith data. 

record (10 char- 

binary record (6 


Content described 

arao ter /word) 

oharacter/word) 


in table 1(b) 

(Cyber code) 

(Field data code) 

2 

First blocked binary 

1000 CDC S.P. word 

1000 Univac D.P. 


record = first 4 

(60 bit) binary 

word (72 bit) 


BET output product 
tape records. Para- 
meter content de- 
sribed in table 1(c) 

record 

binary record 

• 


1000 CDC S.P. word 

1000 Univac D.P. 

• 


(60 bit) binary 

word (72 bit) 

• 


record 

binary record 

N 

Last blocked binary 

1000 CDC S.P. word 

1000 Univac D.P. 


record = up to 4 

(60 bit) binary 

word (72 bit) 


BET output vs. time 
records 

record 

binary record 

EOF 

EOF 


Note: a. Missing quantities will be represented by -9999. 


b. Parameter uncertainty values, except meteorological parameter 
uncertainties will only be included in the accuracy analysis 
output. 
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TABLE I.- Continued 

(b) BET output products tape header record 


OFT number 

Flight launch date: Greenwich year, day (midnight prior to launch) 

Title: 

Ascent Quick-Look BET Output Products 
or Ascent Final BET Output Products 
or Descent Final BET Output Products 
or Ascent BET Accuracy Analysis 
or Descent BET Accuracy Analysis 

Data start and stop time in GMT seconds from midnight prior 
to launch 

Data rate 

Comments 
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TABLE H,- Continued 

(b) Navigation block output tape header record 


78FM*I0 


OFT number 

Flight launch date: Greenwich year, day (midnight prior to launch) 

Title: 

Ascent Quick-Look Navigation Block 
or Ascent Final Navigation Block 
or Descent Final Navigation Block 

Comments 
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TABLE II.- Continued 

(o) Havigatlon block output tape parameter location 


Location 



Number 

of 

elements 

First word 


Last word 

Item 

in 

array 

Record/subscript 

Record/subscript 

Guidance release time 
(4) 

4 

2 

1 

2 

4 

CGT (4,200) 

800 

2 

5 

2 

804 

SET (16,200) 

3200 

2 

805 

6 

4 

REF (27) 

27 

6 

5 

6 

31 

TRAC (6,100) 

600 

6 

32 

6 

631 

MET (20,1000) 

20000 

6 

632 

26 

631 

KAPPA 

1 

26 

632 

26 

632 

PHIO 

1 

26 

633 

26 

633 

DELTA 

1 

26 

634 

26 

634 

LAMO 

1 

26 

635 

26 

635 

RSUBO 

1 

26 

636 

26 

636 

XSUBO 

1 

26 

637 

26 

637 

NAVO(3) 

3 

26 

638 

26 

640 

SPARES 

2360 

26 

641 

28 

1000 
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(d) Navigation 
Parameter 

(lUldanoe release time (4) 
CGT (4,200) 

SET (16,200) 

REF (27) 

TRAC (6,100) 


TABLE II. > Continued 
block output tape parameter description 


Description 


Time when all three IMU's are released (go 
Inertial) , GMT seconds from midnight prior to 
launch 

Shuttle center of gravity history data 


Ti . T2 - T200 

Xi,X 2 X200 

yi, yg ------- - ^200 

Z1.Z2 Z200 


T s time in GET seconds 

X,I,Z, C.G. location (ft) in body axis orthogonal 
coordinates with orgin at nav base (reference 
coordinate A-8) 

Special event tiroes; data will include special 
requested event times 
SET (event, time) 

Event description in 16 words (max. I60 Hollerith 
characters) time of event (max. 200 events) in 
GET seconds (GET from SRB ignition) 

REFSMAT data 

REF1 (3»3) Mean of I960 to IMU-1 stable member 
REF2 (3,3) Mean of 1950 to IMU-2 stable member 
RE1F3 (3,3) Mean of 19&0 to IMU-3 stable member 

Tracker location data, J=1, 100 trackers 
TRAC (1,J) Station ID number 

TRAC (2,J) Station geodetic latitude, deg, + north 
TRAC (3,J) Station geodetic longitude, deg,+ east 
TRAC (4,J) Station altitude, ft, above mean 
sea level 

TRAC (5,J) Station geodetic height, ft, Fischer 
Earth model 

TRAC (6,J) Station azimuth, deg, from true north, 
positive CW 
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TABLE II,- Continued 
(d) Continued 


78PMW 


Parameter 


Description 


MET (20,1000) 


20 meteorological parameters at N points, 


Geodetic latitude, deg, north 
Longitude, deg east to 36 O 
Flag; 0 s measured data, 1 a modeled 
data, 2 s combined measured and 
modeled data 
Spare 

Geometric altitude (above mean sea 
level) in ft 

Horizontal wind speed in ft/sec 
Direction horizontal wind is coming 
from relative to true north, north 
being zero deg, increasing positively 
clockwise in deg 
Ambient temperature in deg R 
Ambient pressure in millibars 
Ambient density in slugs/ft3 
Dew point in deg R 
Ambient temperature systematic 
uncertainty in deg R 
Ambient pressure systematic 
uncertainty in millibars 
Ambient density systematic 
uncertainty in slugs/ft3 
Horizontal wind speed systematic 
uncertainty in ft/sec 
Horizontal wind speed noise or 
fluctuation uncertainty in ft/sec 
Vertical wind speed noise or 
fluctuation uncertainty in ft/sec 
Horizontal wind direction 
systematic uncertainty in deg 
Horizontal wind direction noise 
or fluctuation uncertainty in deg 


N maivimum 

MET 

(1,N) 

MET 

(2,N) 

MET 

(3,N) 

MET 

(4,N) 

MET 

(5,N) 

MET 

(6,N) 

MET 

(7,N) 

MET 

(8,N) 

MET 

(9,N) 

MET 

(10,N) 

MET 

(11, N) 

aMET 

(12,N) 

aMET 

(13, N) 

aMET 

(14, N) 

aMET 

(15, N) 

bMET 

(16, N) 

bMET 

(17, N) 

aMET 

(18, N) 

bMET 

(19, N) 


^Systematic uncertainty is an estimation la error due to instrument errors, 
time and space displacement effects, and any other bias-like factors that 
might contribute to the total low frequency error. 

hNoise or fluctuation uncertainty is an estimation of the la error due 
to high frequency atmospheric inhoraogeneity such as wind turbulence or gusts. 
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TABLE II.- Continued 
(d) Concluded 


Parameter 

Description 

MET (concluded) 

MET (20, N) Spare 

KAPPA 

Plumbline Earth-fixed flight direction from launch 
site (Zg), positive being east of north, deg 

PHIO 

Geodetic latitude of launch site, positive north of 
equatorial plane, deg 

DELTA 

A angle between launch site geocentric radius 
vector and equatorial plane (declination), deg, 
positive north of equator. 

LAMO 

Longitude of launch site, positive east of prime 
meridian, deg 

RSUBO 

Geocentric radius of Earth, at launch site, to the 
vehicle nav base, ft 

XSUBO 

Geodetic height of the vehicle nav base above 
the Fischer ellipsoid, ft 

NAVO 

Structural reference position of nav base (ft) 
in body axis structure-fixed origin 
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Figure 1.- Ascent/descent ancillary support task. 










26 

















78FM40 


appendix 

COORDINATE SYSTQ1 DEFINITIONS AND 
OUTPUT PARAMETER EQUATIONS 


This appendix contains definitions of the coordinate systems listed in table 1(c). 
The definitions were obtained from reference A-1 . In addition to the coordinate 
systen3 definitions, equations for selected table I(o) output parameters are 
given to supplement their tabular descriptions. 


A-1 
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A-1 Davis, Larry D.: Coordinate Systems for the Space Shuttle Program. 

NASA TM X-58153, 1974. 
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NAME: Aries~mean-of-1950 , Cartesian, coordinate system. 

ORIGIN: The center of the Earth. 

ORIENTATION: The epoch is the beginning of Besselian year 1950 or Julian 

ephemeris date 2^133282.423357. 

The plane is the mean Earth's equator of epoch. 

The X|tj axis is directed towards the mean vernal equinox 
of epoch. 


The axis is directed along the Earth's mean rotational 
axis of epoch and is positive north. 

The axis completes a right-handed system. 

CHARACTERISTICS: Inertial, right-handed, Cartesian system. 


Figure A-1.- Aries-mean-of-1950, Cartesian, coordinate system. 


A-3 


I Positfon 

H 



NAME! 

ORIOINi 

ORIENTATION AND DEFINITIONS! 


Velocity 



Plane perpendicular to 


/ ^ 

ArIea-mean“Of-1950i polar, coordinate ayatem. 

For poaibton - the center of the Earth, 

For velocity - the point of intereat, 

P(Xm. VM' Zm5* 

For position, aame as in the Ariea-meBn-cf-1950, 
Cartesian. 

For velocity, reference plane is perpendicular 
ho radius vector Rij from canter of Earth 
to point P of interest. 


Reference direction la northerly along the 
meridian containing P. 


Polar position coordinates of P are as follows, 

right ascension, is the angle between the 
projection of the radius vector in the 
equatorial plane and the vernal equinox of 
epoch, 

declination, is the angle between the 
radius vector and the mean Earth's Equator 
of epoch. 

R|^, magnitude of R).], 

Polar velocity coordinates of P are as follows. 

Let U, E, N denote up, east, and north 
dlreotlons. Then 

azimuth, la the angle from north to 
the projection of on the reference 
plane, positive toward east. 

Y|^, flightpath angle, is the_angle between 
the reference plane and positive 
sense toward U. 

V|^, magnitude of is always positive. 

CHARACTERISTICS: Inertial. 


Figure A-2.- Aries-mean-of-1950, polar, coordinate system. 
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NAME: Greenwich true of date (geographic) coordinate system. 

ORIGIN; The center of the Earth. 

ORIENTATION: The Xq-Yq plane is the Earth's true-of-date equator. 

The Zq axis is directed along the Earth's true-of-date 
rotational axis and is positive north. 

The -*-Xq axis is directed toward the prime meridian. 

The Yq axis completes a right-handed system. 

CHARACTERISTICS; Rotating, right-handed, Cartesian. Velocity vectors 

expressed in this system are relative to a rotating refer- 
ence frame fixed to the Earth, whose rotation rates are 
expressed relative to the Aries-mean-of-1950 system. 

Figure A-3.- Greenwich true of date (geographic). 
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NAME; Csodetio ooordlnnte ayitem. 

OIUGIN: This system consists of s set of paramsters rather than 

a ooordinato syatemi thorefora, no origin is spaolfied. 


ORIENTATION! This system of pararaatero is baaed on an eUlpsoidal model 
of the Earth Ce.g., the Fischer ellipse of i960). For any 
point of interest, we define a Una, known as the geodetic 
local vertical, which <3 perpendicular to the ellipsoid 
and contains the point of Interest . 

h, geodetic altitude, Is the distance from the point of 
intoroat to the rofaronoe ellipsoid, measured along the 
geodetic local vortical, and is positive for points out- 
side the ellipsoid. 

^ Is the longitude measured in the plane of the Earth's 
true equator from the prime (Orecnwioh) meridian to the 
local meridian, measured positive eastward, 

is the geodetic latitude, measured in the plane of the 
local meridian from the Earth's true equator to the geo- 
detic local vertical, measured positive north from the 
equator. 

Note! A detailed explanation of declination, geodetic 
latitude, and geocentric latitude la provided in 
figure A-ll(b). 

CHARACTERISTICS: Rotating polar coordinate parameters. Only position vectors 

are expressed in this coordinate system. Velocity vectors 
should be expressed in the Arles-maan-of-1950, or the Arles 
true-of-date , polar for inertial or quasi-inertial represen- 
tations, respectively, The Fischer ellipsoid model should 
be used with this system. 


(a) Basic definitions. 


Figure A-h,- Geodetic. 
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NAME: Geodetic coordinate system of point P. 

DEFINITIONS: h Is the altitude of point P measured perpendicular from 

the surface of the referenced ellipsoid. 

<t>Q Is the geodetic latitude of point P (angle between the 
true equatorial plane and the geodetic radius of point P). 

4iq Is the geocentric latitude of point P (angle between the 
true equatorial plane and the geocentric radius vector to the 
sub-vehicle point of point P). 

6 is the angle between the geocentric radius vector to point 
P and the equatorial plane (declination). 

X is the longitude of point P angle (•«■ east) between the 
plane of the figure and the plane formed by the Greenwich 
meridian. 

(b) Detailed explanation. 

Figure A-4.- Continued. 
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Let Zbf be the Earth's north polar axis and Xgp pass through the Greenwich 
meridian. Let 


R£ ■ Earth's equatorial radius 

Rp X Earth's polar radius 

e s (Rg>Rp)/R£ X elliptioity or flattening 

To calculate geodetic north latitude, east longitude, Xf and altitude 
above the ellipsoid H, set B a 0.0067 and iterate five times. 

e(2-e)Rg 

° ' ^(*|p * n-e)Zz|p 


Then 


i;. 5 arc tan 


Zef 


V*|f + 'ef 


X s arc tan (Tef/^EF) 


= (' • ® ^ * 2|p 


where Xg^, Yep, and Zgp are Greenwich true-of-date (A-3) values. 


(c) Definitions of 4>* X, H. 
Figure A-4.- Concluded. 
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NAME! Topodstlo coordinate system, 

oniQlN! Orbiter oenter of mass® or navigation base. 

ORIENTATION! is normal to a geodetlo local tangent plane and is 

positive toward the Earth's oenter. 

Xfo is perpendicular to z^d axis and ia positive north- 
ward along the meridian plane containing the Orbiter. 

TTP oompletes the ritct-handed orthogonal system. 

CHARACTERISTICS! Rotating, right-handed, Cartesian system. Velocity vectors 
are expressible in this system for the Oi biter, given 
relative velooity Vto in this system. 


Vtd s (Xjd^ + Yti)3 + Ztd 2)^^2 ^ magnitude of Vtd 


oi; 


•Mi 




VhCr' ' 
05 


Ytd “ sin“^ 



= flightpath angle (relative) 


't’jp' s tan"^ 


»TD 

%D 


= azimuth angle (relative) 


♦d s geodetic latitude 


®A similar system may bo defined for any point of interest. 
Figure A-5.- Topodetlc. 
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The bank aniUi or wigl* of roll about ttM ralatlva wind vaotor la coaputad In 
tha topodatio ayaUa. To daflna iha bank anfla R, P, and V ara flrat 
daflmd raapaotlvaly aa unit vootopa alonn tha roll axta, pitch axla, and wlnd- 
ralatlva valoolty vantor> 

rurthar, daflna U>RxV» (t,a,n)T and a, and n ara dlraoUd ooalnaa of U 
In tha topodatle~ayataa. 



Rotatltn R into V, in tha plana of R and V, daflnaa 0 , and io 
aoooapllshad with aatrlx M (U, 6), alao £' la foraad by tha rotation 



•/ 


■1(R 




M (u, e) 


♦ (1-i2)oo 36) (iB(l-oose)-n slino) (tnO-ooa9)t a alns) 
(ta(1-oon0)+n alns) [a® ♦(1 “b 2 )oos 0) (Bn(1-oos0)-i olnO) 
(tn(1-ooB0)-B aln0) (na(1-oo30)H aln0) (n^ +(1-n2>ooa0) 

£• = M (U,e) £ 


and Vi la a vaotor in tha X-T plana orthogonal to V and la dafinad aa 
Vi a ZXD * V / |ZXD X V| 

finally, bank angla la tha angle between £' and Vi 

Bank Angle = ooa"^ ^£' • VI / |Vi|^ 

and la poaltlve when right wing la below the horizontal. 

Figure A-5a- Concluded. 


A-IO 



NAME: Landing field coordinate system. 

ORIGIN: Runway center at approach threshold. 

ORIENTATION AND DEFINITIONS: Z^p axis is normal to the ellipsoid model 

through the runway centerline at the approach 
threshold and positive toward the center of 
the Earth . 

X^p axis is perpendicular to the Zlp axis 
and lies in a plane containing the Zlp axis 
and the runway centerline (positive in the 
direction of landing). 

Ylf axis completes the right-handed system. 

Alp is the runway azimuth measured in the 
XlfYlf plane from true north to the +Xlf 
axis (positive clockwise). 

CHARACTERISTICS: Rotating, Earth-referenced. 


Figure A-6.- Landing field. 
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/ e raa ^oip2 
[v^h im^sA 

Vvictv' nLssSr^-l^* fromV 

\^S z^V 




^C’ ^c* ^c 

cluster coordinates 
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Note: Xq a X gyro I 

Xq = X gyro I 

Zq a Z gyro I 
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Vertical gyro 


Azimuth gyro 


Rq a Red gyro I.A.J 

Z/^ a Z accel. sens, axis } Single axis accel. 

X^ a X accel. sens, axis \ 

> Dual axis accel. 
a Y accel. sens, axis J 

© a Gimbal torquer 

Gimbal angle resolver 


Figure A-7.- IMU stable member 
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NAME: 

ORIGIN: 

ORIENTATION: 


CHARACTERISTICS: 


Stable Member (IMU) 

The intersection of the innermost gimbal axis and the 
measurement plane of the XY two axis accelerometer. 

The Z(, axis is coincident with the Innermost gimbal axis. 

The Xq axis is determined by the projection of the X 
accelerometer -input axis (lA) onto a plane orthogonal to Zq. 

Yc completes a right-handed triad. 

In a perfect IMU, with all misalinements zero, these 
relationships hold: 

The X accelerometer and X gyro IAS are parallel to the X^ axis. 
The Y accelerometer and Y gyro IAS are parallel to the Y^ axis. 
The Z accelerometer and Z gyro IAS are parallel to the Z^ axis. 

Nonrotating, right-handed, Cartesian system. 

The reference allnement for the gimbal case shall be defined 
with the four gimbal angles at zero and with the vehicle in a 
horizontal position. In a perfect IMU, with all misalinements 
zero and with all gimbal angles at zero, the folic ting 
relationships hold. 

The outer roll axis and the Xq axis will be parallel to Xnb- 
Positive Xq will be in v.he forward direction. Positive roll 
gimbal angles will be in the sense of a right-handed rotation 
of the gimbal case relative to the platform about r.he plus 
outer roll axis. 

The pitch axis and Yq will be parallel to Yjjg. Positive 
Yq will be to the right of an observer looking forward in the 
vehicle. Positive pitch gimbal angles will be in .he sense of 
a right-handed rotation of the gimbal case relati\e to the 
platform about the plus pitch axis. 

The inner roll axis will be parallel to the outer roll axis, 
with the sense of rotation the same as for the outer roll axis. 

The aximuth axis and Zq will be parallel to Zf^g. Positive 
Zq will be down relative to an observer standing in the 
vehicle. Positive azimuth gimbal angles will be in the sense 
of a right-handed rotation of the gimbal case relative to the 
platform about the plus azimuth axis. 

Z|fg are Cartesian components of the navigation base 
coordinate system. 


Figure A-7.- Concluded. 
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ORIGIN; Orbiter center of mass or navigation base. 

ORIENTATION: X^y axis is parallel to the Orbite:' structural body 

Xq axis; positive toward the nose. 

Zby axis is parallel to the Orbiter plane of symmetry 
and is perpendicular to Xby* positive down with respect 
to the Orbiter fuselage. 


Yby axis completes the right-handed orthogonal system. 

CHARACTERISTICS: Rotating, right-handed, Cartesian system. 

L, M, N: Moments about Xby, Yby, aad Zby axes, 

respectively. 

p, q, r: Body rates about Xby, ^BYi and Zby axes, 
respec tively. 

p, q, r: Angular body acceleration about Xby, ^'BY» and 
Zby axes, respectively. 

The Euler sequence that is commonly associated with this 
system is a pitch, yaw, roll sequence, where 6 = pitch, 

'i' = yaw, and 4> = roll. This attitude sequence is pitch, 
yaw, and roll around the Yby, ^by, and Xby axes, 
respectively. 

(a) Basic definition. 


Figure A-8.- Body axes. 
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The wind-relative velocity vector, Veyj, expressed in the body-axis system, 
shall be used bo define the angle of attack, a, and the sideslip angle, 

0, as follows: 



where ^BW, i^BW, ^BW are the components of Vnu in the body -axis system 
(fig. A-8(b)). 


Note: is air-relative velocity; i.e., wind effects are included. 



'BY 


(b) Resolution of wind-relative velocity along vehicle body axes. 


Figure A-8.- Concluded. 
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Earth surface range, S, shall be defined as the product of (1) the magnitude 
of the Earth-fixed position vector to the launch or landing site and (2) the 
central angle between the site vector and the Earth-fixed position vector to 
the Orblter; l.e., 


s * IIsiteI 0E 

RSITE • RoRB 

where Og = COS“i = central angle 

|Rsite| IRorbI 

Rsite = Earth- fixed position vector to launch pad (ascent) or to runway 
coordinate system origin (descent) 

Rqrb = Earth-fixed position vector to Orblter 



Figure A-9 Definition of Earth surface range. 
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Dynamio pressure q shall be defined as 
q = 1/2 p VTd2 

and Maoh number M shall be defined as 

M - VTD 

M - - 5 - 

where VTD = magnitude of wlnd>relatlve velocity 
p = atmospheric density 

a = speed of sound 


Figure A-10.- Definition of dynamic pressure and Mach number. 
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T* 


C s 
00 


(°K) 

Too in -¥ from flight path atmosphere 

V« in fps wind-relative (true) airspeed s VTD 

(N/D) 


/ry 

’t + 122.1 X 10“(5/TJ 
00 

vJ 

T' + 122.1 X 10-<5/T') 

•f 


3.0449939 x 10-8 t 1*5 


M = 


(T + 110.4) 
00 


(lbf-sec/ft 2 ) 




VTDpLB 


(N/D) 


p in slug/ft3 ■+• from flightpath atmosphere 
Lb = body length = 107.5 ft 


= 


VTD VTD 


VTD 


a \4289.05Ta 


(N/D) 


Vi = M 


C„' 


Re 




(N/D) 


Figure A-11.- Hypersonic viscous parameter. 
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Drag over oass (D/M) Is oaloulated as follows: 

, AVs Vrel 

D/M = . r- — r 

At 1 %lI 

where Vrel wind-relative velocity vector 

= ^inertial ” ^inertial 
and b) is the Earth rotation rate 

^inertial position vector 

^inertial = wind-relative Mean of 50 velocity vector 

I'^RElI = magnitude of the wind-relative vector 

AV 3 s sensed change in velocity 

At s time Interval over which AVs is determined 

AVs [y_BD_UNIT X VreJ 
ACC LIFT = . ^ 

At |y bd unit X Vrel I 


Y_BD_UNIT is unit vector passing through Y body (right wing). It is the 
second column of the transformation matrix body to M50. 


LOD = ACC LIFT /(D/M) 


LOD = lift over drag 


LOAD TOTAL = 


AVs 

aT 


/32.I74 


Load total is total load factor. It is total acceleration expressed in g's 
of 32.174 fps2. 

Equivalent airspeed (EAS) = 17.1865 Vi/2 p Vg2 knots 
p = atmospheric density in slug/ft3 
Vg = wind-relative speed in fps 


Figure A- 12.- Calculations for drag acceleration, lift of drag, total 
load factor, and equivalent airspeed. 
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ORIGIN; At the center of the Earth. 


ORIENTATION AND DIRECTIONS: The Xp^ axis is parallel to the gravity gradient 

X' , which passes through the launch site and is 
positive toward the launch site. The XpL axis is 
fixed at SRB ignition. 

The ZpL axis is parallel to and positive in the 
same direction as the chosen Earth-fixed launch 
azimuth direction Z'. 

The XpL axis is parallel to the Y’ and com- 
pletes a standard right-handed system. 

The YpL-ZpL plane is normal to the launch site 
gravity gradient vector. 


CHARACTERISTICS: Inertial, right-handed, Cartesian. 


Figure A-13.- Plumbline. 
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Equations (1), (2), and (3) oaloulata the acceleration due to gravity using the 
eight haraonio gravity equations. 

I 

AGx = -M ^Px ^ ^ - 55 Z - y) in fps2 ( 1 ) 

AGy 5 -M ^Py ^ X - ^ Z - ^ y j In fps2 (2) 

AGz * (pz ^ ^ ^ Z ) in fps2 ( 3 ) 


where 


M 


1.407645W X 10l6 in 


~ 3 el " ® gravitational mass, which Is the 
product of gravitational constant 
and Earth mass. 


Px» Py, Pa, are the total position coordinates in ft 

Initially, the total position value at t-1 s 0 is 

p initial = ECI location of launch pad 

Subsequent first estimates of total position values needed to calculate 
gravitational acceleration is 

/ VSfc + VSt_i \ AGfc.i ^ 

Pt = Pt-1 + y--- — -— + VGt_T y dt + — |-i dt2 in ft (4) 

where 


Figure A-14.~ Gravitational acceleration computation. 
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VS * Mnsfld velocity in fpa 
VG * velocity due to gravity in fps 
(VGt -1 set * 0 at t-1 » 0) 

AG s acceleration due to gravity in fpa^ 
(AGfc.i set * 0 at t-1 s 0) 
dt s tioe Interval in seconds 

P » / Px^ + Py^ ♦ Pz^ in ft 
a * 20925639*8 in ft s radius of Earth 


and 


■ ■ ,x‘ ; t3 (s* • “■ 




(5) 


- ? 2 


8 


Pz 

^ rasi 


(2m- 


DPMMm (c(ram)oosX+a(mm)sinX)m - 


m-t-2 

m+2 


8 


- S —5 [o(nm)cosX+s(nm)sinX ) (-nPn ~ ♦ Pn -1 (n+m)] 
nsffl+1 P“ \ / \ ® /. 


8 

^ - 
m=1 

8 

n=m+1 


(zm2 _ m) ^L— PMM„, ^c(imn)sinX-s(mm)co3X^ 


( 6 ) 


( C(nm)sinX-s(nm)cosX 


\ 


\ ^n*2 

) ^ 


Figure A-l4.- Continued. 
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J 


(Up 
UZ ^ 


Z . I 1211^ FJ(/l) 


fH 

al*Z \cL 


PNPz.i ♦ AXz 


-') 


(7) 


^ Xi_, - (t-1) **.2). 


V ^ m«1 L 


a®*2 / 

p-=^A 


o(an)oosX+s(om)alnX ) 


/ 


ncoH-l 


a*^+2 / p 

n-1) — ( 2 : 

pn+2 \ « 
^oosX c(nm) -f alnX a(nm))j 


(n+m-1) \ 


Where the valuea for the zonal, sectoral and tesaeral harmonic coefficients, 
FJ, and Sp„, respectively, are presented in figure A-16. 


Coefficient selection is accomplished by: 

C(mm)m = C (l+ ^ 10-m\ for m>1, C(mm)i = C(1) 


(8) 


m 


C(nm)„^ * Z C(nm)nm = C(n) (9) 

/ m \ 

S(mm)m s S ll + V 10-m J for m>1, S(mm)i s S(1) (10) 

' i=2 f 

S(nm)nm = 10 -m + n-m^ for m>1, else S(nm)nnj = S(n) (11) 


and 


Figure A-14.- Continued. 
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PNP£ s: ~ PNPa.1 + JKPNji), PNPi - — + 1 


Px^+Py^ 

vm^ s (2m- 1) — P»tln,-1, P^«o = 1 


X = -w»T + tan" 1 (Px/py) in radians = latitude 


Pn = 


Pz 

Pn-1 a (2n-1) - Pn-2 
(n-m) 


Pq = XP, Pi = 0 


Pz Pz 

PNH = PNA-1 — (211-1) - PNa_2<^-1)» PNq = 1» P^l = ~ 


a, 


Pz 

8 ^Pn-1 ot (2n-1) - ^Pn-2 ('^"‘1) 

XP s ^ , XPo = 1, XPi = 

n=1 n 


T = time in seconds from epoch 

0 ) = 0.00007292115147 in radians/sec = rotation rate of Earth 


Figure A-l4.- Concluded. 
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The following equations describe the method of calculating the total position 
and velocity. 

Total position (P) at time t Is defined as the Integral of the total velocity 
from the Initial time to time t. Equation (18) adds the total position 
at time t-1 with the change in position from time t-1 to t. 

/vSt + VSfc_i + VGt + VGt-i\ ^ 

Pt * Pt-1 +l ) 


Where 

p s total position in ft 
VS s sensed velocity In fps (Input data) 

VQ = velocity due to gravity In fps 
Ut = time interval in see 

and 

/ AGt + AGt_i \ 

VGfc = VGt_i + I Jdt in fps (19) 

where 

AG = acceleration due to gravity in fps2 

The velocity due to gravity is calculated by adding the change in velocity 
(due to gravity) to the former velocity (due to gravity). 

Equation (20) describes the total velocity, which is the sum of the sensed 
velocity due to gravity. 

VTt = VSt + VGt ^'Ps (20) 


idiere 


VT = total velocity in fps 

Figure A-15.- Total position and total velocity computation. 
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Zonal Harmon los: 

FJ(I), 1=1,8/ 0.0,0.10826160x10-2, -2.5388x10-6, -1.6560x10-6, 

-2.1848x10-7, 6.2203x10-7, -3.7624x10-7, -3.4040x10-7/ 

Sectoral Harmonica: 


C(I), 1=1,36/ 0.0, 0.0, 2.2181x10-6, -5.0409x10-7, -8.8795x10-8, 

-5.3098x10-8, 2.5437x10-7, -6.8551x10-1°, 1.5765x10-6, 3.1196x10-7, 
7.6894x10-8, 1.0551x10-7, 7-6476x10-9. 3.1555x10-8. 5.5912x10-9, 
9.8324x10-8 5.8692x10-8, -1.4740x10-8, 1.0411x10-9, 3.5092x10-9, 

1.1576x10-1°. -4.0641x10-9, -2.1795x10-9, -2.8471x10-1°, 
-5.8546x10-1°, -3.2546x10-1°, 4.^52x10-1°, -2.0509x10-'°, 
6.7546x10-12, -5.7460x10"12 1.7718x10-12 -2.4693x10-11, 

-1.7303x10-12, -3.1326x10-13, 3.3836x10-13, -1 .6308x10“l3/ 

Tesseral Harmonics: 

S(I), 1=1,36/ 0.0, 0.0, 2.8843x10-7, -4.5586x10-7, -1.0960x10-7, 

2.3455x10-8, 1.0557x10-7, 3.9182x10-8,-9.0602x10-7, -2.2055x10-7, 
1.4562x10-7, -5.1423x10-8, -4.3734x10-8, 1.1015x10-8, 3.8638x10-9, 
1.9611x10-7, -1.1854x10-8. -6.9874x10-9, 7.7530x10-1', 

-2,9593x10“9, -7.6414x10-1°, 6.7006x10-9, 4.6325x10-1°, 

-1.7714x10-9 -2.6051x10-10, 7.3081x10-11, -1.6231x10-9, 
-4.4558x10“1°, 1.1795x10-11, 1.4325x10-11, -5.6889x10-1', 
1.0979x10-11, 8.5920x10-12, 2.4477x10-13, 4.3704x10-13, 
1.5645x10-13/ 


Figure A-16.- Harmonic coefficients. 
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